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Применение математической регрессии для определения параметров 
трёхфазного асинхронного двигателя” 
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(Донской государственный технический университет) 


Статистический разброс параметров электротехнических устройств часто является проблемой для анализа 
их работы и точного определения оптимальных режимов. В этих УСЛОВИЯХ требуются математические алго- 
ритмы, позволяющие производить диагностику таких устройст, В С ТОЧНОСТЬЮ, достаточной для инженерных 
расчётов. В настоящей работе предложен новый алг Оритм расчёта характеристики момент-скольжение 
трёхфазног о асинхронного двигателя. Он позволяет определить величины максимального момента и крити- 
ческого скольжения. Для этого статистически обрабатывается начальный участок характеристики, получен- 
НЫЙ экспериментальным путём. Этот в значительной степени универсальный метод решения подобных элек- 
тротехнических задач базируется на общих математических подходах. В связи с этим он может быть поло- 
жен в основу электронного диагностического уст| ройст, ва на базе микропроцессора. 

Ключевые слова: математическая регрессия, метод наименьших квадратов, характеристика момент- 
СКОЛЬЖЕНИЕ, диагностика электротехнических уст. ройст, В. 


Введение. Применение сложной микропроцессорной техники позволяет создавать устройства 
диагностики электродвигателей, основываясь на минимальном количестве измерений. Их нехват- 
ка компенсируется сложными математизированными алгоритмами. На сегодняшний день это не 
повышает стоимости подобных приборов, но значительно упрощает для потребителя процедуру 
диагностики. Таким образом, соответствующее направление представляется актуальным. В 
настоящее время сформировались некоторые теоретические и практические подходы к вопросам 
диагностики электродвигателей, в том числе асинхронных [1—3]. Одно из таких решений пред- 
ставлено в настоящей статье. 

1. Постановка задачи. Как известно, основной энергетической характеристикой трёхфазного 
асинхронного двигателя является характеристика момент-скольжение (рис. 1), выражаемая с до- 
статочной точностью формулой Клосса [4]. 

м (5) = я (1) 

кр 

Как видно из (1), для её получения необходимо знать значения двух величин: максимального мо- 
мента М, и критического скольжения $5„,. Эти величины известны не всегда. Их эксперимен- 


тальное обнаружение затрудняет остановка двигателя при достижении им точки [5 Мы)- При 


теоретическом нахождении указанных величин сложность представляет определение параметров 


Т-образной схемы замещения. Таким образом, выяснив значения величин М, и 5,„,, можно с до- 


статочной точностью рассчитать все остальные точки рабочего участка кривой М5). 
2. Описание метода. Предлагаемый в данной статье метод основан на статистической обработ- 
ке начального участка характеристики М5) по методу наименьших квадратов (МНК) [5] с вычис- 
лением величин Ми $,„,. Как следствие, получаем теоретическую зависимость (1). 





* Работа выполнена в рамках инициативной НИР. 
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Как известно [5], МНК основан на определении экстремума некоего функционала 
я . 
Р(алаьниа,) = У Уже, — Улей (Хит )| >тт, (2) 


где у, — вектор экспериментальных значений функции, у "| — вектор теоретических 

(предполагаемых) значений этой функции с аргументами, полученными экспериментальным 
путём. 
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Рис. 1. Нормированная характеристика момент-скольжение 




































































Сделаем некоторые преобразования в формуле (1): 
2 2 
1 9 +59 _ 5 бе 
М (5) 2М„5.5 2М„5„ 2М„5 





Введём обозначения: аа 5, =Ь, - =У(5). 

















М, ы М (5) 
Тогда 
а5 ‚ ав 
: 3 
7) 55 +55 (3) 
Из (2) следует, что 
2 
а5, ар | 
пеьь)- Уи, а. > тм. (4) 
Тогда в точке экстремума 
0! (а,Ь) аз, а 5 В 
= | 1 1 =0, 
да > 2Ь = 2Ь р (5) 
ов а, аб | 25 а |_| 
05 9) Ь 25 |2 25“ 
После преобразований получим 
1 и 
[5 И : 7.5 | ау. -0 ей 





[БЕ а-я иь-о 
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Окончательно для определения неизвестных коэффициентов а и Ь получим следующую систему 
нелинейных уравнений: 


а аа 
ребре + ана + аа" =0, 


Ь 
Ва а? (7) 
пря + Вьа + оз +Вна?Ь =0, 


где коэффициенты ал и В„ определяются согласно (6), п — количество измерений. 

Эта система может быть решена численно, например в приложении Ма САО, что и было 
проделано в настоящей работе. 

С практической точки зрения интересно изучить поведение полученного алгоритма в 
условиях вносимых погрешностей измерений при определении величин скольжения и вращающе- 
го момента. 

В этой связи после успешного тестирования алгоритма при точных значениях величин $, 


и М, эти величины брались с отклонением в 3 % от их истинных значений. Испытания проводи- 


лись с различным количеством измерений (от 7 до 50), а также с различными интервалами по 
шкале величин скольжений ($ < 0,55, $ <0,65,, $5< 0,75). 


кр! кр! 
На рис. 2, 3 представлены результаты этих исследований для величин математических 


ожиданий М и среднеквадратичных отклонений о ошибок определения искомых величин 5р и 


М, . Величины 5, брались с одинаковыми интервалами. 





























































































































М, % 5 < 0,55 
5 
4 5х 0,65 

+» о — — 
й т | ьь > 1 ыы В Е =_= 
\ 

1 
0 5 < 0,75 

0 10 20 30 40 50 “ 


Рис. 2. Зависимости математического ожидания (в процентах) ошибки измерения от количества измерений на различных 


интервалах по шкале величин скольжений 
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Рис. 3. Зависимости среднеквадратичного отклонения (в процентах) ошибки измерения от количества измерений 


на различных интервалах по шкале величин скольжений 


87 


Технические науки 








Данные рис. 2, 3 иллюстрируют рост устойчивости в определении искомых величин при 
увеличении объёма выборки и ширины участка по шкале величин скольжений. Начиная с объёма 
выборки, равного 40, устанавливаются удовлетворительные (не превышающие значения вноси- 
мой погрешности) значения величин Ми о. Однако для кривой, соответствующей $ < 0,55,„,, эти 


значения недопустимо велики, поэтому целесообразно применять интервалы измерений 
5 <0,65„, ИЛИ $ <0,75„. 

Поскольку устойчивость расчётов определяется в том числе величиной интервала по оси 
скольжений, то представляется желательным получить возможность определения этого интерва- 


ла каким-либо алгоритмическим способом. 
На рис. 4 представлены нормированные (к величине максимального момента) зависимости 


т(5), а также её первая и вторая производные 71’(5) и т"(5) соответственно. Видно, что функ- 
ция 77"($) имеет минимум в интересующей нас области, который можно использовать для опре- 


деления необходимого интервала по оси скольжений. 


4т(5)/ 46?'103 





Рис. 4. Зависимости нормированных функций момент-скольжение, а также её первой и второй производной 
от скольжения 
М (5) 
М 


т 





Для нормированной функции момент-скольжение т(5)= перейдём к новой величине 


5 
5, =— , тогда 
кр 











т (5, ) = РА, (8) 
$+— 
1 
Трижды дифференцируя (8), получим: 
(5) [1 Е ВВ И ВЕ 
1 1 ] 52 1 ] —_ 52 1 ] я : 
$5+— 5$+— $+— | $, 
1 5. 5. 
Решая уравнение л7” (5, } = 0, получим: 
4 _ 2 
12. 51 65, 80. 
$2 +1 
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Отсюда 
5-65? +1=0. 
Это уравнение имеет корни: 
5; = (-2,414 -0,414 0,414 2,414), 


из которых условию задачи удовлетворяет $, = 0, 414 , то есть 5„„ = 0,4145. 


Значит, отследив равенство нулю третьей производной функции л7(5), можно с некото- 
рой точностью «масштабировать» ось скольжений, определив таким образом необходимый для 
измерения интервал. Учитывая, что погрешность численного дифференцирования растёт с ростом 
степени производной, в данном случае лучше следить за функцией л7”(5), определив величину 


5 Как аргумент, при котором эта функция имеет минимальное значение. 


ти 
3. Результаты и выводы. Данный метод может быть положен в основу диагностического 
устройства, работающего на базе микропроцессора. Для диагностики можно использовать и ПК, 
считывающий информацию непосредственно с измерительных датчиков. При этом, если отсут- 


ствуют какие-либо данные о величинах $, и М, то для их поиска исслед емый двигатель дол- 
у кр т! 


жен сделать два прогона. Первый — ДЛЯ приближённого определения величины интервала по оси 
скольжений. Для этого нужно найти $, ЧТО ПОЗВОЛИТ ОПределить величину шага $, . На втором 


прогоне будут найдены параметры $5„, и М,. 


Для реализации этого метода использовались общие математические подходы. Поэтому 
сфера его применения не исчерпывается решением представленной здесь задачи. Так, например, 
модифицируя регрессионную функцию, с его помощью можно экспериментально определить па- 
раметры однофазного трансформатора Р, и Р. , задавшись известной в электротехнике формулой 
для зависимости КПД от коэффициента загрузки [4]: 


_ В$,, со$ф 
п) В.5, с05ф +Р, + В?Р ` 
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